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Abstract. 

The time-dependent Dirac equation is solved using the three-dimensional Finite 
Difference-Time Domain (FDTD) method. The dynamics of the electron wave packet 
in a vector potential is studied in the arrangements associated with the Aharonov- 
Bohm effect. The solution of the Dirac equation showed a change in the velocity of 
the electron wave packet even in a region where no fields of the unperturbed solenoid 
acted on the electron. The solution of the Dirac equation qualitatively agreed with the 
prediction of classical dynamics under the assumption that the dynamics was defined 
by the conservation of generalized or canonical momentum of the electron. 
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1. Introduction 

The Finite Difference- Time Domain (FDTD) method, originally introduced by Kane 
Yee p] to solve Maxwell's equations, is for the first time applied to solve the time- 
dependent three-dimensional Dirac equation. The Zitterbewegung and the dynamics 
of well-localized electron were used as examples of FDTD applied to the case of free 
electrons [2]. The motion of electron wave packets inside and scattering from the 
potential step barrier or linearly dependent potential, arrangements associated with the 
Klein paradox [3] , were used as examples of electron interaction with an electromagnetic 
scalar potential pfj. In this paper, a FDTD study of the dynamics of an electron wave 
packet under the influence of a vector potential is presented. Such a dynamic behavior 
is most often associated with the Aharonov-Bohm effect [SJ. 

In the Aharonov-Bohm effect, as a manifestation of quantum mechanics, charged 
particles passing around a long solenoid can feel a magnetic flux even when all the 
fields of the unperturbed solenoid are zero in the region through which the particles 
travel. The shifts in the phase of the wave functions describing the particles have been 
experimentally verified by its effect on the interference fringes [HI El El E] • Since for the 
unperturbed solenoid, there are no classical forces acting on the charged particles in 
the zero field region, the theoretical description of the Aharonov-Bohm effect contains 
a number of assumptions. They include assumptions on nonlocal features of quantum 
mechanics, the physical meaning of the vector potential, topological effects, etc., but 
generally accepted physical understanding is still lacking [JU]. While there is still an 
open question on the presence of classical forces responsible for the Aharonov-Bohm 
effect PIU IU H2], on a macroscopic level they have not been observed [T3] . 

Proper quantum-mechanical description of the dynamics of a relativistic charged 
particle involves the solution of the Dirac equation in the time domain. In addition 
to initial conditions, such a dynamics is defined only by the configuration of the 
electromagnetic scalar and vector potentials, and does not involve knowledge or any 
assumption on "classical forces". The solutions of the time-dependent Dirac equation, 
some of which are presented in this paper, can shed light and fill critical knowledge gaps 
on the theoretical and experimental interpretations of the mechanism of the Aharonov- 
Bohm effect, including the existence or non-existence of classical forces. 

2. Time-dependent solution of the Dirac equation 

The FDTD solutions of the time-dependent Dirac equation were obtained for the case 
when the electromagnetic field described by the four-potential = {A (x), A(x)} was 
minimally coupled to the particle [TU [15] 
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The matrices a and (5 were expressed using 2x2 Pauli matrices a s and the 2x2 unit 
matrix I. 

In the FDTD method, the time dependent solution of the Dirac equation was 
obtained using updating difference equations [2]. As an example, the values of at 
the position (iAx,jAy, kAz) and at the time step (n+ 1/2) At were obtained using the 
equation 
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where C n (I, J, if) = 1 + ij^[mc 2 + eA^I, J, if)]. The space and time were discretized 
using uniform rectangular lattices of size Ax, Ay and Az, and uniform time increment 
At. While it is not generally required, in the Eq. (|5| Ax = Ay = Az. Updating 
equations for \l/ 2 , \l/ 3 , and ^4 were constructed in a similar way. as a result, the dynamics 
of a Dirac electron can be studied in any environment described by a four-potential A^ 
regardless of its complexity and time dependency. In this paper we studied the dynamics 
defined by the vector potential A(x) ^ 0. 

The Dirac equation is a differential equation of the first order and linear in d/dt. 
As in the case of Maxwell's equations, the entire dynamics of the electron is defined, 
only by its initial wave function. In this paper, the dynamics of a wave packet was 
defined by the initial wave function of the form 
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where N = [(27r) 3/,2 Xo] _1//2 is a normalizing constant. Eq. (|6j represents a wave packet 
whose initial probability distribution is of a normalized Gaussian shape. Its size is 
defined by the constant xo, its spin is pointed along the z-axis, and its motion is defined 
by the values of the momenta p\, p%, and p^. Some consequences of the initial localization 
of the wave packet on the overall dynamics of the electron were studied in Ref . [21 US] ■ 
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3. Validation of the computation 

3.1. Dynamics of a wave packet in a strong uniform magnetic field 

The dynamics of a wave packet is very complex. The dynamics of a particle described 
by the wave packet in Eq. ^ depends on its localization, defined by the Gaussian 
component of the wave function, and its initial momentum, part of the wave function's 
phase. While an extensive study was done on the dynamics of the wave packet related 
to the scalar component A (x) of a four-potential [4], such a study does not validate the 
dynamics of the wave packet related to the vector component of a four-potential, A(x). 

Applying the FDTD method to study the dynamics of a wave packet in a vector 
potential associated with a strong uniform magnetic field is not difficult. Classically, the 
dynamics consist of uniform rotational motion. In the relativistic quantum-mechanical 
description, however, even such a simple dynamics can validate the FDTD computation 
only up to a certain level. As pointed out in Ref. [17J "the dynamics is particularly 
rich and not adequately described by semiclassical approximations" . In Ref. [18] it was 
demonstrated that in the presence of an external magnetic field the wave packet splits 
into two parts which rotate with different cyclotron frequencies, and after a few periods, 
the motion acquires irregular character. As a result of such a complex dynamics, when 
comparing the results obtained by the FDTD method and the computation described in 
Ref. [T7j and |18j . we could not expect more than a qualitative agreement. In addition 
to this qualitative agreement, of equal importance to the validation of the computation 
is the consistency of the results obtained for different vector potential gauges. 

In this paper, the motion of the wave packet described by Eq. ^ was studied in a 
uniform magnetic field oriented along the y-axis 



For this field the corresponding vector potential in a rotationally invariant gauge is 



The dynamics of the wave packet in a uniform magnetic field was first obtained by 
solving the Dirac equation for this vector potential. 

While the probability densities \^>\ 2 were calculated for the entire computational 
volume and at every time step, their values are shown here only in the horizontal plane 
or plane of classical particle motion. The schematics of the wave packet or charged 
particle motion relative to the orientation of the magnetic field is shown in Figure [TJ 
The position and the shape of the wave packet in the horizontal plane, as it moves along 
its first orbit, is shown in Figure [2] The initial position of the wave packet was at the 
center of the vector potential. Its initial momentum was p± = 0.53 MeV/c, making the 
motion relativistic. In order to force a relativistic electron to complete a full circle in the 
available computational space, the magnitude of the magnetic field was Bq = 10 s T. The 
field of such a strength is associated with the fields at the surface of the neutron stars. 
The classical orbit of the electron in this field is r c i ass = pi/(eB Q ) = 1.76 x 10~ 2 nm. 



B = (0,fl ,0). 



(7) 




(8) 
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Figure 1. Schematics of the classical orbit of the electron relative to the orientation 
of the uniform magnetic field. 

As shown in Figure |2j during the first rotation the wave packet generally follows the 
classical orbit, but disperses at the same time. The position of the center of probability 
of the wave packet relative to the classical orbit is shown in Figure [3] 

The complexity of the dynamics of the wave packet increased in the later stages of 
the motion. While the wave packet followed the circular motion, the probability density 
|^| 2 at some times assumed the spiral shape shown in Figure 4 increased or reduced its 
length, changed its rotational motion, and translated from one place to another. Overall, 
the characteristics of this dynamics were similar to the characteristics described in Ref. 
[TT] and [18]. While the figures show some of the complex shapes of the probability 
density |\I/| 2 , the richness of the electron motion can be better appreciated through the 
animation of the dynamics accessible on-line |19j . 

Several additional tests of the FDTD method were performed using the dynamics 
of the electron in the uniform magnetic field. As expected, due to the normal spin and 
the magnetic field orientation, reversal of the orientation of the magnetic field resulted 
in the reversal of the direction of the rotation of the wave packet, keeping the same 
properties of the dynamics of motion. 

Of particular interest was testing the effects of the choice of gauge. The motion of 
the same wave packet was also studied for the uniform magnetic field oriented along the 
y-axis defined by the translationally invariant gauges 

A = B Q (z,0,0), (9) 



or 

A = B (0, 0, —x). (10) 




Figure 2. The shapes and the positions of the wave packet in the horizontal plane 
during the first rotation. The circle represents the classical orbit. The animation can 
be accessed on-line [T9] . 
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Figure 3. Positions of the center of probability of the wave packet in the horizontal 
plane during the first rotation. The circle represents the classical orbit. The animation 
of the dynamics can be accessed on-line [19]. 



In both cases the wave packet persisted in a circular motion following the classical 
orbit. While the dynamics of the wave packet behaved as expected, using translationally 
invariant gauges has an additional importance on validating the FDTD computation. 
As seen in Eq. (|5]), the gauge in Eq. ^ couples to the ^ 4 component and the gauge 
in Eq. (10) couples to the ^ 3 component. Similarly, the cross coupling exists for other 
components of the wave function As a result, the same dynamics should be obtained 
by different combinatorics of the components of the vector potential and the components 
of the wave function. This enabled for testing of possible inconsistencies in the FDTD 
updating equations. No inconsistencies were found. 

The dynamics of the wave packet motion in a uniform magnetic field was used here 
only as a validation of the FDTD method when a vector potential was applied in the 
Dirac equation. The same complexity of the quantum dynamics of motion as in previous 
publications p21 HE] was shown. While used here only for computational validation of 
the FDTD method, this dynamics could be studied as a separate problem in more detail 
in the future. 

Finally, the complexity of the quantum dynamics of the wave packet in a uniform 
magnetic field studied here for three choices of gauge could be fully appreciated only by 
downloading the related animations [T9] . 
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Figure 4. Changes in the shape of the probability density |^| 2 during one of the 
subsequent orbits of the wave packet. The circle represents the classical orbit. The 
animation of the dynamics can be accessed on-line |19j . 

3.2. Motion of a wave packet between two infinite lines of electric dipoles 

The motion of a wave packet in a uniform magnetic field is used to validate the FDTD 
method when a vector potential A(x) is applied in the Dirac equation. The motion of a 
wave packet between two infinite lines of electric dipoles is used to validate the dynamics 
of a wave packet under influence of a position-dependent electromagnetic scalar potential 
Aq(x), and to validate the geometry associated with the Aharonov-Bohm effect studied 
in the next section. In the case of a wave packet moving between two infinite lines of 
electric dipoles, the dynamics obtained from the solution of the Dirac equation can be 
compared to the dynamics of a charged particle influenced by a classical electromagnetic 
force. The classical dynamics, resulting from the interaction of a charged particle and 
two infinite lines of electric dipoles, is obtained following Boyer [20J. 

The scalar potential A Q of two infinite lines of electric dipoles parallel to the y-axis 
and separated by a distance 2a, with dipoles oriented in the x-direction, parallel or 
anti-parallel to the direction of the wave packet motion, is 
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Figure 5. On the left side, the solution of the Dirac equation with the scalar potential 
Aq described by Eq. (11) is compared to the solution of Eq. (12 1. The right side 
compares the solution of the Dirac equation with the scalar potential described by Eq. 
(13) to the solution of Eq. (14). The red circles rapresent the solutions of the Dirac 
equation and the solid lines the classical solutions. In both cases the momentum of 
the charged particle was p — 0.09 MeV/c at the position x — —0.4 nm. The distance 
between the lines of electric dipoles was 2a = 0.38 nm and the line density of dipoles 
was p = 1.7 10~ 17 Cm/m. In the coordinate system used in this paper two infinite 
lines of electric dipoles were parallel to the y-axis and the wave packet moved along 
the x-axis. 



where p is the line density of dipoles and eo is the permittivity constant. With the choice 
of z = 0, the classical motion of a charged particle along the x-axis, under influence of 
this potential, can be calculated by solving the equation 

dp x _ ^dA _^ qp 1 / 2x 2 \ ^ 
dt dx Txe^x 2 + a? \ x 2 + a 2 J 

The comparison of the quantum-mechanical dynamics obtained as a solution of 
the Dirac equation with the scalar potential described by Eq. (11) and the classical 
dynamics obtained as a solution of Eq. ( 12 ) is shown on the left side in Figure [5] The 
relativistic effects in the classical dynamics were reduced to less than 2% by choosing 
the initial momentum p = 0.096 MeV/c. The lines of electric dipoles were separated by 
a distance 2a = 0.38 nm and the line density of dipoles p was 1.7 10~ 17 Cm/m. Figure 
[5] shows agreement between the quantum-mechanical and classical dynamics over the 
distance for which the solution of the Dirac equation was computed. 

Under the same conditions, the computation was repeated for dipoles oriented in 
the z-direction, anti-parallel to each other, and perpendicular to the direction of the 
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wave packet motion. In this case the scalar potential Aq is 

a =+ P ( z - a z + a \ no) 

2vre V^ 2 + «) 2 x 2 + {z + a) 2 J ' 1 ' 

and the equation of motion, for z = 0, is 
dPx = 2ax 

rft T 7re (a; 2 + a 2 ) 2 ' 1 ' 

As shown on the right side in Figure [5j the solution of the Dirac equation with the scalar 



potential described by Eq. (13) and the solution of Eq. (14) agree again over the range 
of computation. 

4. Wave packet dynamics in a vector potential created by two infinite 
solenoids 

The goal of this work is to study the dynamics of an electron wave packet under 
the influence of a vector potential associated with the Aharonov-Bohm effect [5]. 
Particularly, in this paper we present the dynamics of a wave packet obtained from the 
solution of the Dirac equation with a vector potential created by two infinite solenoids. 

The vector potential of a single infinite solenoid oriented along the y-axis can be 
written as 



[z, 0, —x) for r < R 



a = < ' ' ' - " , (15) 

to?s{z,0,-x) forr>R 

where $ = Bq-rRq is the magnetic flux inside the solenoid, B defines the strength of 
the magnetic field, Rq is the radius of the solenoid, and r = \/x 2 + z 2 is the distance 
from the center of the solenoid in the x-z plane. Outside two parallel infinite solenoids 
separated by a distance 2a and with opposite magnetic field orientation, the vector 
potential is then 

a ^ ( z—a z+a x x \ 

= 2^\x 2 + {z- a) 2 ~ x 2 + (z + a) 2 ' ' x 2 + (z + a) 2 ~ x 2 + (z - a) 2 ) ' ( ' 

An example of the shape of this potential is shown in Figure [6j Outside the solenoids, 
the associated magnetic fields B = V x A of the vector potentials described by Eqs. 



(15) and (16) are zero. 

The time-dependent Dirac equation is now solved for the scalar component of the 
four-potential A = and the vector component A defined by Eq. jl6| . With this 
choice we can study the dynamics of the wave packet in the region where the electric 
and the magnetic fields of the unperturbed solenoids are zero. 

The dynamics of the electron motion between two infinite and parallel solenoids is 
essentially the same as the electron dynamics in the case of Aharonov-Bohm effect. In 
both cases the electron moves in a region where the unperturbed solenoids produce no 
Lorentz force on the electron. Because of the symmetry, however, the dynamics in the 
case of two infinite solenoids consists only of the motion along one straight line between 
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Figure 6. The shape of the vector potential of two parallel infinite solenoids with 
oppositely directed magnetic fluxes of the same magnitude. The potential is shown in 
the plane normal to the orientation of the solenoids. In the coordinate system used in 
this paper the magnetic field of the solenoids is oriented along the y-axis, the vector 
potential lays in the x-z plane, while the wave packet moves along the x-axis. 

the solenoids. This avoids the complications of the Aharonov-Bohm dynamics where the 
trajectories on opposite sides of a single solenoid are compared. Because the properties 
of the dynamics are not studied through an interference, in the case of two solenoids 
the changes in the dynamics of the electron wave packet exclude the contribution of the 
nonlocal properties of the electron wave function dependent on the topology of the space. 
Another advantage of using two solenoids and localized wave packet is the possibility 
to separate the solenoids far enough that no parts of the wave packet penetrate the 
non-zero field region inside the solenoids, excluding this contribution to the dynamics as 
well. If necessary, the solenoids could be shielded by a potential barrier of a supercritical 
potential additionally preventing any penetration of the wave packet into non-zero field 
region. 

The motion of the wave packet between the solenoids, in the plane normal to 
the orientation of the solenoids, is shown in Figure [7j The wave packet was initially 
positioned away from center of the solenoids. The dynamics was studied for two initial 
momenta, p\ = 0.53 MeV/c or p\ = 0.64 MeV/c. The solenoids were separated by a 
distance of 2a = 0.1 nm and the magnitude of the magnetic flux inside each of them 
was $ = 5.2 x lO" 14 Wb. 

As shown in Figure [7], and in the related animation, the wave packet moved between 
two solenoids along a straight line, dispersing in time in the direction normal to the 
direction of motion. Since for the unperturbed solenoid, there was no classical force 
acting on the electron, one should have expected a constant velocity of the wave packet 
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Figure 7. Six stages of the motion of the wave packet between two infinite parallel 
solenoids. Small yellow circles represent the positions and the sizes of the solenoids. 
The animation can be accessed on-line [T§] , 



along the straight line. This was not the case. The velocity of the wave packet, shown in 
Figure [8] as a function of the position, increased as the packet approached the solenoids 
and decreased as the packet left the solenoids. It is somewhat paradoxical that while 
there was no classical force expected to act on the electron, the electron acceleration, 
obtained strictly as a solution of the Dirac equation, was not zero. 

Extensive studies and verifications of this result were performed. They include 
the dependence of the solution of the Dirac equation on particle charge, energy, spin 
orientation, orientation of the solenoids, and the strength and the orientation of the 
magnetic field inside the solenoid. In addition, the time and the position dependencies 
of the particle momentum, energy and canonical momentum were computed. The results 
could be summarized as follows: 

• the change of the velocity of the wave packet did not depend on the particle spin 
orientation relative to the orientation of the solenoids 

• as shown in Figure [9j change of the sign of a charged particle inverted the order of 
the change of the velocity, but the velocities converged to the same value at spatial 
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Figure 8. The velocity of a wave packet normalized to the speed of light as a function 
of the wave packet position. The wave packet is moving along a straight line between 
two infinite parallel solenoids positioned at x=0. 



infinity 

• change of the orientation of the magnetic field inside the solenoids had the same 
effect as change of the sign of charged particles 

• as also shown in Figure [9j the mechanical momentum of the particle changed during 
the process, but the total energy and the canonical momentum did not change 



5. The dynamics of the charged particle under the conditions of the 
Aharonov-Bohm effect 

The change of velocity was predicted more than 35 years ago by Liebowitz [2T] and 
by Boyer [22], and was proposed as the basis for the Aharonov-Bohm phase shift. It 
was shown that if the magnetic interaction energy of the unperturbed solenoid currents 
and the passing charge is assigned to a kinetic energy change of the passing electron, 
then the electron has the same qualitative velocity change as found from our numerical 
calculation of a wave packet dynamics in the Dirac equation. Through lowest order 
in the particle-solenoid interaction, the energy conservation criterion used by Boyer is 
equivalent to conservation of the passing charge's canonical momentum p c = p + qA. 
Thus 

= ±[(p^ + m 2 cY 2 + qv-A] 

pc 2 dp _ -> -> 
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Figure 9. Left: change of the velocity of a wave packet normalized to the speed of 
light as a function of the wave packet position. Top curve rapresents electron, bottom 
curve positron. The red dots depict the positions of the solenoids at x=0. Right: 
the top curve shows the change in the electron wave packet momentum. Mid and 
bottom lines show that the total wave packet energy and canonical momentum are 
conserved. The wave packets were moving along a straight line between two infinite 
parallel solenoids separated by 0.072 nm. The magnitude of the magnetic flux inside 
each of them was $ = 4.2 x 10~ 14 Wb. Initial momentum of the wave packets was 0.53 
Mev/c. 



_ d(p + qA) _ dp c 

dt dt K ' 

where the term q(dv/dt) ■ A was neglected as higher order in the interaction. This 
equation implies conservation of the canonical momentum of the particle. In classical 
dynamics, the motion of an electron can now be obtained by solving the differential 
equation 

d(p + qA) 



dt 



0. 



Since A x does not explicitly depend on time, for the motion along the x-axis, Eq. (18) 
can be written as 

dp x _ _ ( dA^ _ _^dA^(ix^ _ _ ( dA^ v 
dt dt dx dt dx x 
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Figure 10. Higher, red curve, represents the velocity, as a function of the position, 
of the relativistic particle moving along a straight line between two infinite parallel 
solenoids, obtained as the solution of Eq. (22 1. Lower, black curve, represents the 



velocity of a wave packet obtained from the solution of the Dirac equation under the 
same conditions. The velocities are normalized to the speed of light. The solenoids are 
positioned at x=0. 



In the relativistic case 

dp x _ 
dt 



m 



dv x 
dt 



(20) 



Here mo is the rest mass of the electron. Using Eqs. (18) and (20) we get the differential 
equation for the velocity of the electron 



q 



dv x 

dt mo 
Substituting A x from the Eq. 



dA x 
dx 



(21) 



(16) for z = 0, electron motion along the mid-path 



between two solenoids, the differential equation for the velocity of the electron becomes 



dv n 



2g$ 



ax 



dt 7rmo (x 2 + a 2 ) 2 ^ c 2 J ' ^ 
where a is half the distance between two parallel infinite solenoids. This differential 
equation can be solved numerically. Figure [lO] shows the qualitative agreement of the 



classical dynamics of the particle obtained from the solution of Eq. (22 ) and the quantum 
dynamics of the corresponding wave packet obtained from the solution of the Dirac 
equation. 

To conclude this section, we can look at the ongoing question if the Aharonov-Bohm 
effect is the result of a force changing the velocity of a particle passing on opposite sides 
of a single infinitely long solenoid, or there is only a quantum-mechanical phase shift 
[231 EU [25], |26l E7J EH]. Quantum-mechanically such a question does not exist. The 



Dynamics of the Aharonov-Bohm effect 



16 



dynamics of the relativistic electron should be obtained only as a solution of the time- 
dependent Dirac equation. The solution of the Dirac equation shows that the velocity 
of a wave packet changes even in the region where the unperturbed solenoid gives zero 
magnetic field. Since the change of the velocity depends on the gradient of the vector 
potential, the velocity of the wave packet passing on the opposite sides of the solenoid 
will be different. However, the particle exits the near-solenoid region with the same 
energy as when it entered, regardless of the side on which it passed the solenoid. 

From the studies presented in this paper, the classical picture of the Aharonov- 
Bohm effect consists of the phase shift in the wave function of the particles passing on 
opposite sides of a solenoid being attributed to time lag resulting from different evolution 
of the velocities of the particles. A phase shift then results in an interference pattern 
change. 



6. Conclusion 



In conclusion, the full three-dimensional Finite Difference Time Domain (FDTD) 
method was developed to solve the Dirac equation. In this paper, the method was 
applied to the dynamics of the electron wave packet in a vector potential in an 
arrangement associated with the Aharonov-Bohm effect. The solution of the Dirac 
equation showed that the velocity of the electron wave packet changed even in the 
region where the electric and the magnetic fields were zero. 

The solution of the Dirac equation qualitatively agreed with the prediction of 
classical dynamics under the assumption that the dynamics were defined by the 
conservation of generalized or canonical momentum. 

The studies in this paper help establishing a picture of the Aharonov-Bohm effect 
as the interference pattern resulting from the phase shift in the wave function of the 
particles passing on opposite sides of a solenoid attributed to a time lag resulting from 
different evolution of the velocities of the particles. 
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